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Abstract. - We study heat conduction in one dimensional (ID) anharmonic lattices ana- 
lytically and numerically by using an eftective phonon theory. It is found that every eftective 
phonon mode oscillates quasi-periodically. By weighting the power spectrum of the total heat 
flux in the Debye formula, we obtain a unified formalism that can explain anomalous heat con- 
duction in momentum conserved lattices without on-site potential and normal heat conduction 
in lattices with on-site potential. Our results agree very well with numerical ones for existing 
models such as the Fermi-Pasta-Ulam model, the Frenkel-Kontorova model and the 4>'^ model 
etc. 



Recent years has witnessed increasing studies on heat conduction in one dimensional (ID) 
anharmonic (nonlinear) lattices [1]. On the one hand, people would like to know whether 
or not the Fourier's law of heat conduction for bulk material is still valid in ID systems. 
This is a fundamental question in non-equilibrium statistical mechanics. In fact, it is not 
trivial at all as a rigorous proof is still not possible. On the other hand, the fast development 
of nano technology makes it possible to fabricate ID or quasi ID systems such as nanowire 
and/or nanotube etc and to measure its transport properties. To understand heat conduction 
behavior in such systems is of great interest in heat control and management at nanoscale. 
Numerically, an anomalous heat conduction - heat conductivity diverges with system size- 
has been observed in momentum conserved systems without on-site potential such as the 
Fermi-Pasta-Ulam (FPU) lattice [2], and a normal heat conduction has been found in the 
systems with on-site potential like the Frenkel-Kontorova model [3] and the 0"* model [4,5]. 
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Unfortunately, up to now a general theory to predict the heat conduction behavior in a ID 
system is still lacking. 

In this Letter, we investigate analytically and numerically the physical mechanism leading 
to the anomalous and the normal heat conduction in ID anharmonic (nonlinear) lattices from 
an effective phonon theory. The theory is based on the ergodic hypothesis (equipartition 
theorem). As it will be seen that, our analytical result can explain the anomalous and normal 
heat conduction observed numerically in different models. 

We consider a ID anharmonic (nonlinear) lattice with the Hamiltonian 



N 



i=l 



rriixl + V{5xi^i+\) + U{xi) 



(1) 



with periodic boundary condition xi = xn+i- Here we chose mass mi = 1 for all lattice. 
Xi = dxi/dt. 6xi^i^i = Xi — Xi^i. Without loss of generality, we write the inter particle 
potential V{6xi^i^i) and the on-site potential U{xi) as 



V{6xi,i+,)=J2i 



{Sxi^i+iY 



s=2 



s=2 



(2) 



respectively. The canonical transformation which diagonalizes the harmonic Hamiltonian is 
X = BQ, where Xi = {xi oiXi), Qk = [qu oipk), and B^k are [6] 



Bik = 



tGfcCos^^^ fc = l,...,[f] +1 
|G,sinHM^™) k=[f]+2,...,N 



where [^] is the integer part of ^, and Gk = l/V^ for = 1 and k = N/2 + 1 if iV is even, 
otherwise Gk = 1. The spectrum of harmonic lattice is Uk =2 sin {k — 1)it/N. 

Under the ergodic hypothesis, the system obeys the generalized equipartition theorem [7] 

fcsT = ^9fef^^, here (•) denotes the canonical ensemble average. The force in k space has 
two parts 



dH 



N oo 



ffe = ^ = X! X] {'^kgsiSxi^i+iY ^^ik + a^Xi ^Bik) , 



= 1 s=2 



where the new matrix 7jfe is defined as in Ref. [6], 7jfe = for fc = 1 and 7,^ = {l/LOk){Bik — 
Bi+ik) otherwise, ^'.k satisfy X^ili Sxi^i+ilik = ^kQk, Yl,k=2lik'^kqk = The general- 

ized equipartition theorem becomes 



N 



ksT = ^^{(jJkgsiiSxi^i+iY ^qk)lik+(^s 



[xt ^qk)B^k) 



=1 s=2 



E 



9s ,^^r ,^ —'^k + '^s 



J2^=li^Xi,i+lY) 



O^i^k + l){(lk) 



(ql) 



(3) 



where 



J2T=2 9AY.f=iiSxi^i+iy) 



a 



{Elii5xi,i+,r} 
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and 

In analogy with an Harmonic lattice where ksT — uff.{q^), we define effective phonons in 
ID anharmonic lattices. The frequencies of the effective phonon are 

(4) 

here k is wave vector with replacement k - 
velocities are 

"'^^^ 

It has been found [6, 8] numerically that the approximation in Eq. Q is feasible for 
anharmonic lattices without on-site potential such as the FPU-/? model {U{x) — 0, V{x) = 
jl -X- fix"^ we take /? = 1 in this paper) and model iJJix) = 0, V{x) — a;'*/4.), a is 
found to be independent of mode k and lattice length. 

In order to check if Eq. Q is suitable for the anharmonic chains with on-site potential, 
we study numerically 7 for the (^"^ model {U{x) — x'^/2,V{x) = x*/A) and the quartic cf)'^ 
model {U{x) — x^/A, V{x) = x'*/4). We find that 7 is independent of k and the chain length 
provided the chain is long enough. For example, 7 ~ 1.065 for the (j>* model with chain length 

> 64 at temperature T = 1. In Fig^ we show {Pk) / i'^kilk)) ^^"^ dispersion relation dik for 
quartic (j)'^ (upper panel) and (j)'^ (lower panel) models. {p1) / i'^li^k)) '^^^y close to one for 
all k in both models. We also find that 7 increases with temperature as 7 '--^ j^o.eiio.oi ^ggg^ 
Fig. 121). Therefore, we can conclude that the approximation is good for both anharmonic 
chains with and without on-site potential under ergodic hypothesis. 

It should be pointed out that the approximation in Eq. is a kind of mean-field approx- 
imation. It is based on our numerical observation as a rigorous proof is not possible,. The 
phonon-phonon interaction is implicitly contained in the two coefficients a and 7. 

The heat conductivity can be derived from the Debye formula 

k 

where Ck, Vk, and are specific heat, phonon velocity and phonon relaxation time of mode 
fc, respectively. Generally, the contribution from phonons of different frequencies should in 
principle be weighted but this is not reflected in the Debye formula. Since the heat conductivity 
is a transport coefficient of heat energy which manifest itself by heat flux, this weight factor 
must be the power spectrum of total heat flux. Thus we rewrite the Debye formula for ID 
anharmonic lattices as, 

K = ^ PfcCfeWfcTfe, < A: < 27r, 

k 

in which is the normalized power spectrum of total heat flux. Moreover, the Debye formula 
does not give an explicit expression of phonon relaxation time Tk . From our numerical calcu- 
lations, we find that each mode oscillates quasi-periodically for both lattices with and without 
on-site potential (see FigO. It is therefore reasonable to assume that the phonon relaxation 
time is proportional to the quasi-period of each mode, i.e., r^, = A3^, where prefactor A is 
only temperature dependent and will be discussed in another paper. 



27r(fc — 1)/N and ujk = 2 sin |. The corresponding 



cos-. (5) 
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Fig. 1 Fig. 2 

Fig. 1 - The energy ratio, {Pk) /^>k{lk) ^ f'lid dispersion relation uik versus k for the model (lower 
panel) and the quartic (f)'^ model (upper panel). A'^ = 128 for the (f)'^ model, and A'^ = 64 for the 
quartic model. Temperature is set to T = 1. 

Fig. 2 - Parameter 7 versus temperature T for the 0'' model. 



To deal with the size dependence of the heat conductivity, it is more convenient to consider 
the mean-free-path rather than the relaxation time. The mean-free-path of the effective 
phonons is defined by 

Ik = VkTk = 2nX—^ — cos - (6) 

^fc + 7 2 

For systems without on-site potential where 7 = 0, reduces to 

27rA k 
Ik = cos-. 

uJk 2 

In the long wave-length limit, fc — > 0, the mean-free-path cx 1/fc, becomes divergent. 
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However, for systems with on-site potential where 7 > 0, the mean-free-path of any effective 
phonon is finite. This difference resuhs in different heat conduction behaviors in systems 
without on-site potential and systems with on-site potential. 

The modified Debye formula of thermal conductivity can be expressed in a continuous 
form in the thermodynamical limit: 

K = — P{k)ckVklkdk = cXy/^ / P{k) cos^ -dk (7) 

Jo Jo (ujI+j)^ ^ 

where c = X^fc^fc. By definition, P{k)dk — P{uj)duj, and P{lu) is the Fourier transform of 
auto-correlation function of total heat flux J{t). 

Eq. (UJis the main analytical result of this Letter. Whether the heat conduction in a ID 
system is normal or anomalous depends on whether the integral of Eq. Q is finite or infinite. 
For systems with on-site potential such as the FK model and the model where 7 > 0, 
all integrands except the normalized power spectrum P{k) are finite. Since P{k) should be 
normalized over the phonon spectrum J^^ P{k)dk = 1, the integral of Eq. {Tj) for systems 
with on-site potential is always finite. Thus even without knowing exact knowledge of P{k), 
we can predict that heat conduction of systems with on-site potential obeys Fourier's Law. 
However for systems without on-site potential, 7 = which implicitly means the momentum 
conservation. The integral of Eq. {Tj) reduces to J^^ P{k) ''"^^ ^ dk. In the long wave-length 
limit, ojk ~ k, the integral has a singularity at fc ^ 0. This singularity originates from the 
infinite mean- free-path of effective phonon of the long wave-length limit and is proportional to 
1 /k. Since the effective phonon of long wave-length is the most dominant part for heat transfer, 
the systems without on-site potential will exhibit an anomalous heat conduction. In the 
following we shall apply above theory to two typical classes of Hamiltonian systems, the FPU 
model, a representative model without on-site potential, and the (f)^ model, a representative 
one with on-site potential. As shall be seen soon that our theory gives predictions in good 
agreements with numerical simulations. 

Systems without On-Site Potential. In the Hamiltonian U{x) —0 means momen- 
tum conservation, and 7 = in Eq.©. Here we focus on the FPU-/? model. The details of 
renormalized frequencies of effective phonons will be discussed as well as the effective phonon 
speed and the heat conductivity with respect to lattice length. 

In the FPU-/3 model, V{x) = {2 + x^/^ and U{x) = 0. a is determined only by 
temperature [6,9]. In this model, a has a simple analytic expression 

this equation is equivalent to the Eq. 1 1 of Ref . [9] . 

If the effective phonons of long wave-length are dominant in heat transfer, which is true 
for a system of large size, the speed of energy transport can be approximated by the speed of 
effective phonon with the longest wave-length 

-y==^/acos^. (9) 



In Fig. 21 we draw above analytical v versus temperature for two different system sizes and 
compare them with numerical ones from Aoki and Kuznezov [10]. The agreement between our 
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Fig. 4 



Fig. 3 - Time behaviors of qk{t) of the H4 model and the quartic 
both models, = 2it /uk- 



model. iV = 64 and T = 1 for 



Fig. 4 - The effective phonon speed in FPU-/3 model vs temperature. The numerical results come 
from the lower panel of Fig. 2 in Ref. [10] 



analytical results © and numerical ones are very good in a very wide range of temperature 
( five orders of magnitudes.) 

In the thermodynamic limit N 00, the size dependence of heat conductivity only de- 
pends on the integral 



2tt 



P{k) 



9 h 

COS^ -I 

^dk. 



(10) 



/o t^fc 

If the effective phonons of long wave-length dominate the heat transfer, and P(fc) cx 
{6 > 0) asymptotically in long wave-length limit, then, the heat conductivity is a system size 
dependent quantity, k cx oc . Numerical calculation of power spectrum of total heat 
flux in FPU-/3 model has shown that S « 0.37 ~ 0.4 [2]. This divergent behavior of thermal 
conductivity has been observed in the FPU model by different groups [2,4, 10]. 

Systems with On-Site Potential Things turn out to be different for systems with on- 
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site potential. Having an on-site potential means that 7 > and the momentum conservation 
is broken. A phonon band gap appears around zero frequency. Under the ergodic assumption, 
heat conduction in such systems obeys the Fourier's law. For the sake of simplicity, we only 
discuss the cj)^ model. 

The 4>'^ model has a quadratic inter particle potential, V{x) — jl and a quartic external 
potential U{x) = Therefore, we have a = 1 and 7 = 2;f)/(X]i=i 2^?)- The 

ensemble average has no simple analytic expression. At temperature T = 1, the numerical 
value 7 « 1.065 is found to be independent of length at least for N > 64. The dispersion 

relation, Cj^ — \J 4 sin^ ^ + 1.065, is shown in Figure ^ 

The mean-free-path depends on temperature via parameter 7 and A, Ik = v^t^ = 
2A7r sin fc/(4 sin^ | -t- 7). In the low temperature limit, the (j)^ model reduces to the harmonic 
model where A is infinity. Therefore, A should decrease with temperature. From Fig. 12 we 
can see that 7 increases monotonically with temperature. As a result, the mean-free-path 
decreases as temperature is increased. Since the heat conductivity is independent of N when 
N is larger than the mean-free-path. The higher the temperature, the smaller the size effect 
for the (f)^ model as observed by Aoki and Kuznczov [10]. 

In summary, from an effective theory, we have derived an analytic formula for heat con- 
ductivity in ID nonlinear lattices, Eq. lO. We find that the phonon-phonon interaction due 
to the anharmonicity (nonlinearity) can be written as an effective harmonic one in terms of 
ensemble average, so we can attribute the heat transfer to the effective phonons which can be 
treated in the same way as phonons in the Harmonic lattice. The difference between system 
without on-site potential and system with on-site potential lies in their renormalized frequen- 
cies. For systems without on-site potential, there exists zero frequency mode with infinite 
mean-free-path which is the physical mechanism for anomalous heat conduction. For systems 
with on-site potential, a phonon band gap near zero frequency appears leading to the normal 
heat conduction, i.e. heat conduction obeys the Fourier's Law. 
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